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Abstract 

Astrophysical observations and cosmological data have led to the conclusion that nearly one quarter of the Universe 
consists of dark matter. Under certain assumptions, an observable signature of dark matter is the annual modulation of 
the rate of dark matter-nucleon interactions taking place in an Earth-bound experiment. To search for this effect, we 
introduce the concept for a new dark matter experiment using Nal scintillation detectors deployed deep in the South 
Pole ice. This experiment complements dark matter search efforts in the Northern Hemisphere and will investigate 
the observed annual modulation in the DAMA/LIBRA and DAMA/Nal experiments. The unique location will permit 
the study of background effects correlated with seasonal variations and the surrounding environment. This paper 
describes the experimental concept and explores the sensitivity of a 250 kg Nal experiment at the South Pole. 
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1. Introduction 

An abundance of astrophysical and cosmological 
observations have established a modern concordance 
model of cosmology, known as ACDM 1 1 1. This model 
prescribes that 23% of the Universe can be attributed to 
cold dark matter, 4% to baryonic matter and the remain- 
ing 73% to dark energy. Much of what we know about 
dark matter is inferred solely from its gravitational inter- 
actions. While evidence for dark matter has been firmly 
established |2|, dark matter has yet evaded direct obser- 
vation and thus many of its aspects remain speculative. 

Weakly interacting massive particles (WIMPs) are a 
theoretically favored candidate for dark matter |3 |. A 
suite of direct detection experiments is now underway 
worldwide to search for WIMPs through observation 
of WIMP-nucleon elastic scattering. The sensitivity of 
these experiments is now sufficient to allow interesting 
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theoretical regimes to be probed, such as those predicted 
by the Minimal Supersymmetric Model [4]. 

A WIMP moving at typical galactic rotation speeds 
of ~220 km/s could in principle interact coherently with 
a nucleus, inducing a nuclear recoil with energy of a 
few tens of keV |i5J. Detecting such interactions re- 
quires stringent control of background and long, sta- 
ble exposure times in low-noise environments. This 
can be achieved by operating detectors in ultra-low- 
radioactivity environments and at underground loca- 
tions to minimize activity induced by cosmic rays. 
Many such detectors make use of phenomena by which 
the particle interaction energy is observed through more 
than one channel (e.g. phonon and ionization), the ratio 
of which allows nuclear recoil events to be distinguished 
from an electron recoil background |6 |. 

An additional strategy is to search for a modulation 
in the detected event rate. Under commonly accepted 
models, galactic dark matter is distributed in a diffuse, 
roughly spherical halo through which the solar system 
travels as the Sun moves within the disk of the Milky 
Way|7J. The rate of events detected by an Earth-bound 
experiment depends on the relative velocity of the dark 
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matter to the detecting medium. The rotation of the 
Earth about the Sun produces a small but periodic varia- 
tion in this relative velocity. This affects the rate of dark 
matter interactions in the earth-bound detector and the 
energy spectrum of the induced recoils. If the coupling 
to baryonic matter is large enough, a periodic change in 
the rate of particle interactions may be observable lH. 
Such a modulation would exhibit a period of one Earth- 
year and peak at the time of year when the Earth-dark 
matter relative velocity is at an extremum. 

The DAMA/Nal experiment and its successor, 
DAMA/LIBRA, have been operating an array of sodium 
iodide (Nal) crystals in the Laboratori Nazionali del 
Gran Sasso (LNGS). They have reported a statistically 
significant observation of an annual modulation, which 
has been attributed to dark matter ||9] [TOl . However, the 
interpretation of these data as evidence for elastic scat- 
tering of WIMPs off sodium (Na) or iodine (I) nuclei re- 
mains inconsistent with the observations from a number 
of direct detection experiments IfTTHTTll . Moreover, the 
presence of a signal below 6 keV, as claimed by DAMA, 
would require a significant suppression of background 
continuum compared to energies above 6 keV, in con- 
tradiction with Monte Carlo simulations of radioactive 
background ifTSl . 

The challenge of reconciling the DAMA observations 
with those of other experiments has generated a great 
deal of interest in alternative models of dark matter 
(e.g. |T9'-'2T|). To date, none of the proposed scenarios 
has satisfactorily accounted for the DAMA modulation. 
It remains that the DAMA claim has never been tested 
at an independent location with the same detector target 
material (Nal) and an experimental setup with equiva- 
lent sensitivity. Given these motivations, we describe 
an annual modulation search using Nal crystals that can 
be deployed in the South Pole ice. 

2. Annual modulation search at the South Pole 

2.1. Conceptual Description 

This paper describes a Nal dark matter experiment 
that may be deployed to a depth of 2450 m in the Antarc- 
tic icecap. This Southern Hemisphere location pro- 
vides a unique opportunity to disentangle seasonal ef- 
fects from an annual modulation caused by dark mat- 
ter. The ice offers a quiet and stable operating environ- 
ment, a necessary requirement for an annual modula- 
tion search. The proposed location is at the bottom of 



suggested location meets these goals. Details of back- 
ground estimates and sensitivity are provided in subse- 
quent sections. 

An annual modulation arising from dark matter will 
have the same phase regardless of the location on Earth. 
Seasonal effects, however, are opposite in the Northern 
and Southern Hemispheres. For instance, at the South 
Pole, the peak of a dark matter signal that is consis- 
tent with the DAMA result will occur six months out 
of phase from the peak in cosmogenic muon flux and 
associated background. Such a location thus allows one 
to disentangle seasonal environmental effects from as- 
trophysical effects. Further, the depth of ice provides 
a significant overburden that is comparable to under- 
ground laboratories, significantly reducing the flux of 
cosmogenic muons passing through or near the dark 
matter detectors. Data from IceCube provide detailed 
measurements of the muon flux. Correlations with these 
data may be exploited to remove any remaining cosmo- 
genic background from a potential dark matter signal 
(Section|23]i. 

An experiment consisting of 250 kg of ultra-pure 
Nal crystal scintillators, which is comparable in size 
to DAMA/LIBRA, can confirm or exclude the DAMA 
dark matter claim with a minimum of two years oper- 
ation time (Section |3]l. Key to such an effort is the de- 
velopment of crystals that have a comparable or lower 
internal radioactive contamination than the DAMA Nal 
crystals. The DAMA/LIBRA experiment observed a 
single-scatter event rate of ~1 count per day (cpd) 
per kg per keVee (electron-recoil equivalent energy) in 
the 2-lOkeVee region of interest |T0|. The observed 
modulation amounts to a ~2% change in the rate of 
events observed between 2-6keVee- Reproducing such 
a measurement requires achievement of a very low en- 
ergy threshold and good control of environmental back- 
ground. 

The Antarctic ice provides a radiopure and stable en- 
vironment for a dark matter experiment. Radioactivity 
levels in Antarctic ice are estimated to be comparable 
to the cleanest shielding materials used in underground 
experiments, removing the need for additional shield- 



the IceCube neutrino detector. Section 2.1 provides an 
overview of the basic design and requirements for such 
an experiment. It also describes the way in which the 



ing from the ice (Section 2.2 1. The level of radon in the 
ice is expected to contribute a negligible component of 
the total background. The temperature of the ice at the 
South Pole varies in time by less than a degree. Based 
on a known profile measured by IceCube, the tempera- 
ture at the surface is -50° C and increases to -20° C at a 
depth of 2500 m. The measured and calculated profiles 
of the temperature versus depth at various sites across 
Antarctica agree to within 0.2° [22 1. At such tempera- 
tures, the light output of a Nal crystal is several percent 
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higher compared to room temperature; however, scintil- 
lation will occur over a longer time frame Il23ll24l . 

With the successful installation and operation of 
nearly 5200 digital optical modules in ~ 1 km-^ of ice, 
IceCube has demonstrated the technical feat of safely 
deploying delicate instrumentation in Antarctic ice and 
has achieved reliable, remote operation of these de- 
vices. The scientific infrastructure at the Amundsen- 
Scott South Pole Station, operated by the National Sci- 
ence Foundation, provides excellent technical and sci- 
entific support. Utilizing much of the technology and 
expertise developed by IceCube, it is possible to at- 
tach encapsulated Nal crystals and their PMT systems 
to deployment-readout cables similar to those used in 
the IceCube array. A variety of deployment depths, up 
to a maximum of 2500m, are under consideration. 

Challenges associated with a deep-ice deployment, 
although not insurmountable, must be adequately ad- 
dressed in the experimental design. Dark matter detec- 
tors require a method of robust in-situ calibration. Pos- 
sibilities within the ice include LED's, and natural and 
cosmogenic radioactivity in Nal. Transport to and stor- 
age at the South Pole Station, at an altitude of 2850 m, 
will result in cosmic activation of the Nal crystals, read- 
out devices and encapsulation assembly. Depending on 
the activated isotopes, this activity will decay away on 
a variety of timescales and may result in a period of 
higher background at the start of the data collection pe- 
riod. Once a detector is deployed, it becomes perma- 
nently inaccessible and thus reliability of operation will 
be critical. Nal scintillators represent an ideal choice 
for such a requirement because they are relatively sim- 
ple to operate and have a history of use in field-based 
research. These issues and others are currently under 
study for optimization in a polar-ice deployment. 

2.2. Background from Environmental Radioactivity 

Gammas, betas, alphas and neutrons will arise from 
trace radioactive contaminants in the crystals and the 
surrounding material (PMTs, pressure vessel, shielding, 
ice, etc.). Some of these particles will produce signals 
in the Nal that are indistinguishable from nuclear recoils 
expected from WIMPs. The annual modulation search 
is a compelling strategy for detectors that cannot dis- 
tinguish between nuclear and electron recoils. In such 
cases, one may search for a modulation in event rate. 
As discussed in Section [3] the sensitivity of such an ex- 
periment is limited by the level of background that is 
achieved. 

We perform an initial estimate of such background 
events expected in a l-lOkeVge window. The estimate 



consists of a Geant4 simulation ll25l of radioactive de- 
cay from several primary components of a detector as- 
sembly and the surrounding ice. We consider a unit con- 
sisting of a single 8 kg crystal that is instrumented with 
two light guides, two PMTs and encased in a stainless 
steel pressure vessel (PV) as shown in Figure [T] Typi- 
cal levels of radioactivity for the assembly materials are 
input to the simulation and are summarized in Table [T| 




Figure 1 : Geant4 visualization of tlie simple detector assembly used 
for radioactive background estimates. Shown is the steel PV (large 
cylinder with enlarged cap on the left side and rounded bottom on the 
right side), with an 8 kg Nal crystal inside (small, interior cylinder). 
The crystal is instrumented with a light guide and PMT on each end. 



Table 1: Assumed concentrations of U, Th and '""K, in ppb for 
major components of a Nal assembly and the surrounding ice. Details 
on the estimate for contamination in Antactic ice are in the text. 



Material 








drill ice ||23 


0.076±0.046 


0.47±0.14 


<262 


Antarctic ice 


10-4 


10-4 


0.1 


PMT 1261 


30 


30 


60000 


steel PV II23 


0.2 


1.6 


442 


Nal 


0.005 


0.005 


10 



In the simulation, the ice surrounding the detector is 
treated as two distinct volumes, the drill ice and the 
glacial ice. In reality, the volume immediately sur- 
rounding the pressure vessel is expected to consist of 
a ~60 cm diameter column of "drill" ice that is formed 
when hot, drill water mixes with melted glacial ice and 
refreezes in the drill hole. Surrounding this hole will be 
pure, Antarctic glacial ice. To assess the activity levels 
expected from the drill ice, samples from several holes 
were taken during IceCube construction. Using 2 liter 
polypropylene bottles, the samples were gathered from 
the drill return water, which was pumped out from the 
top of the drill hole. The samples were counted at SNO- 
LAB fTT\. The results of the count are shown in Ta- 
ble [T] Improvement of the contamination level can be 
achieved by changing the recirculation of water through 
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Figure 2: The histograms show the simulated activity in Nal crys- 
tals due to uniform radioactivity in the ice surrounding the detector 
assembly. 



the drill. We simulate the drill ice contamination within 
a column of ice with a diameter 40 cm larger than the 
diameter of the vessel and a length 2 m longer than the 
vessel. The results are shown in Table |2] 

The glacial ice background is simulated by assuming 
the contaminants are uniformly distributed in a spher- 
ical volume of ice, 2 m in radius, surrounding the de- 
tector. Figure |2] shows the simulated activity in the Nal 
crystals arising from contamination in this volume of 
ice at the level of 1 ppt ^^**U and ^^^Th and 1 ppb of 
""'K. We scale the activity levels shown in Figure [2] to 
match the contamination levels expected in the Antarc- 
tic glacial ice. Details on the estimated radiopurity of 
the Antarctic ice are in the next paragraph. The simula- 
tion results are shown in table|2] We note that this calcu- 
lation overestimates the glacial ice background because 
we have already accounted for the volume immediately 
surrounding the pressure vessel with the drill ice sim- 
ulation; however, the contribution of the glacial ice is 
nearly negligible so this assumption has little impact on 
the results. 

Although no ice core samples have been extracted at 
the South Pole, we can estimate the glacial ice contam- 
ination levels from the measurements taken at Vostok 
located 1300 m away from the South Pole lH ||29]| . 
Most of the contaminants arise from volcanic ash ac- 
cumulated as dust layers in the ice. Optical measure- 
ments from AMANDA and IceCube allow us to corre- 
late the dust layers present at the South Pole to those 



at Vostok 1221 iSOj. The ice at the South Pole between 
2100m to 2500 m deep is roughly 70,000 to 100,000 
years old, and contains about 0.1 ppm of dust. Vol- 
canic dust typically contains about 1 ppm of ^^'^U and 
^^^Th and lOOOppm of ""'K [311, yielding contamina- 
tion levels of ~10 "^ppb in the ice for ^^^U and ^^^Th 
and 0.1 ppb for ""'K. The ^^^U concentration was mea- 
sured independently in the Vostok ice core samples us- 
ing inductively coupled plasma sector field mass spec- 
trometry (ICP-SFMS), yielding consistent results ESll . 

The results of the simulation for all considered mate- 
rials is summarized in Table |2] Assuming a radiopurity 
for Nal crystals that is comparable to DAMA, the net 
predicted rate of events from the considered sources is 
1.3-1.7 cpd/kg/keVee in the 1-lOkeVee window where a 
WIMP signal is expected. Although not shown in the 
tables, we also considered ^"^Co contamination of the 
steel vessel at the level of 7.2 mBq/kg |27|, which adds 
an additional 0.1-0.2cpd/kg/keVee to the total. The es- 
timated background due to radioactivity is thus compa- 
rable to the total event rate observed by the DAMA col- 
laboration. Despite the assumption of high radiopurity, 
the Nal crystals themselves remain a large source of ra- 
dioactivity. Care must be taken to minimize this contri- 
bution as much as possible. In this hypothetical setup, 
the steel vessel (Sandvik SANMAC SAF 2205 hollow 
bar duplex steel) is a significant component and can be 
further reduced by constructing the encapsulation from 
cleaner material or implementing a lining within the 
pressure vessel to shield the crystals. 

Table 2: Shown are the estimated contribution to event rate from 1- 
lOkeVec in a single 8 kg Nal crystal. The first three items are cal- 
culated using the Geant4 simulation of a simple 8 kg assembly. The 
internal Nal contamination is taken from ( 1 8 1 where the energy spec- 
trum of events from radioactivity was simulated for the DAMA exper- 
iment assuming radioisotope concentrations reported in |32|. 

Material event rate in Nal 

(cpd/kg/keVee) 

drill ice 0.8 

Antarctic ice < 0.001 

photomultiplier tubes 0.01-0.02 

steel PV 0.2-0.6 

Nal crystal ~0.3 



Neutrons, which will produce nuclear recoils in the 
Nal crystals, are also generated as a result of radioac- 
tivity through fission and (a,n) processes. In these es- 
timates, we neglect this source of background because 
the rate of neutrons generated by radioactivity is several 
orders of magnitude lower than the rate of gamma-rays. 
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Such events will thus be a small component of the net 
rate observed in the Nal scintillator 



being designed, will enhance the effectiveness by sev- 
eral orders of magnitude. 



2.3. Background from Cosmic Rays 

Cosmic muons produce background for dark matter 
experiments in the form of spallation products. These 
consist of neutrons and excited nuclei, which may be 
long or short-lived. Recently, long-lived phosphores- 
cent states have also been postulated to be a product 
of cosmic activity in Nal crystals |33|. Owing to the 
large overburden achieved at underground laboratories, 
background produced by cosmic muons is expected to 
be well below the rate of events caused by radioactiv- 
ity. However, the cosmic muon flux, and thus the activ- 
ity induced by it, is modulated with a period of about 
one year due to the change of temperature and pressure 
in the upper atmosphere. Seasonal variations will thus 
play an important role in modulating this portion of the 
event rate and must be thoroughly understood. 

It has been widely noted that the maximum of the 
DAMA signal (May 26 ±7 days QOl) and the peak of the 
muon flux at LNGS (July 5±15 days f34]) bear a close 
proximity in time. The amplitude of the muon modu- 
lation at the South Pole is 10% and is strongly corre- 
lated with the atmospheric temperature. The maximum 
muon rate occurs in mid-January [35J. This modulation 
is larger than the 2% eff'ect at LNGS (EH. Thus, any 
unremovable eff'ects of cosmic ray origin that can be at- 
tributed to the DAMA signal will be opposite in phase 
and more pronounced at the South Pole. Correlations 
between minor fluctuations in temperature and muon 
rate have also been observed and are well-characterized 
with IceCube data ||35l . Such a high degree of correla- 
tion offers the opportunity to reject a cosmic-ray origin 
for any potential modulating signal. 

It is worth noting that, at a depth of 2450 m in the ice, 
the rate of muons and any associated spallation prod- 
ucts will be comparable to the rate at deep underground 
sites. Using MUSUN Eg), the muon intensity at 2200 
m.w.e. is estimated to be 1.3 x lO^-'yu/m^/s withamean 
energy of 307 GeV/c^. Ice is more eff'ective than rock 
at moderating low-energy (<10MeV) neutrons down to 
energies where they are unable to produce a nuclear re- 
coil above threshold in a Nal detector. Additionally, any 
remaining short-lived activity that happens to propagate 
into a crystal will occur in coincidence with a muon and 
thus can be vetoed. IceCube has already demonstrated 
its utility as a robust atmospheric muon veto. Within 
the volume of DeepCore, a low-energy extension to Ice- 
Cube, a reduction of 8 x 10"^ in muon rate has been 
achieved ||37l . Advanced software algorithms, currently 



3. Sensitivity 

The DAM A/LIBRA and DAMA/Nal experiments are 
currently the world's most sensitive Nal-based dark 
matter experiments. It is expected that an array of 
Nal detectors with the same mass, energy threshold and 
intrinsic background level as DAMA will achieve the 
same level of sensitivity as the DAMA experiments; 
however, it is less clear what level of sensitivity can 
be achieved with detectors that have a different reach 
in background level and energy threshold. To explore 
these parameters, we calculate the sensitivity of a series 
of simulated experiments. 

We consider an annual modulation arising from spin- 
independent elastic scattering of WIMPs, which are dis- 
tributed according to the standard halo model, by fol- 
lowing the prescription in Il38l . Under these assump- 
tions, the modulation of events at a specific recoil en- 
ergy can be approximated with a cosine function as fol- 
lows: 



dR(t, E) 
dE 



So(E) + S,„(E)cos(oj(t-t,)), 



(1) 



where oj - Itt/T and T = I year. The value tc is the 
time of year when the velocity of the WIMP halo rel- 
ative to the Earth-bound experiment is at a maximum. 
The energy, E, corresponds to the target recoil energy. 
The sign of the modulation, S„i, can be positive or neg- 
ative and depends on kinematic constraints. S[) is the 
unmodulated component of the dark matter signal. 

The DAMA collaboration analyze their data by mea- 
suring a residual rate of modulating events over the ob- 
served constant rate in a fixed interval of energy. We 
follow the same procedure and integrate over the energy 
interval which corresponds to the region of experimen- 
tally accessible recoil energies: 



> Q(m) 



(E)dE 



^max ^thresh 

This yields a simplified modulation equation: 

R(t) - Sq + S,„cos(aj(t - tc)). 



(2) 



(3) 



For the case of multiple target species in one detector, 
as with Nal, the net recoil rate is the sum of the recoil 
rates off' the individual nuclei, weighted by the respec- 
tive mass fractions. Because the time dependence re- 
mains the same for different targets, the simplified form 
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Figure 3: The curves show the sensitivity of hypothetical 500kg-year exposures with varying total event rates (in cpd/kg/keVec). Shown are 
two energy threshold scenarios. The left plot shows sensitivities with a 2 keVee experimental threshold. The right plot shows sensitivity with a 
4keVec threshold. The gray regions show the 90% (dark) and 99.7% (light) DAMA/LIBRA allowed regions for interactions with Na (masses of 
~ 10 GeV/c^) and I (masses of ~ 100 GeV/c^). DAMA/LIBRA allowed regions are calculated without channeling. 



for the time dependent rate shown in equation[3]remains 
unchanged. 

Experimental data are typically displayed using the 
calibrated electron-recoil equivalent energy (keVee) be- 
cause Nal cannot provide event-by-event discrimina- 
tion between nuclear and electron recoils. Working in 
the same units, we present results in two intervals in 
electron-recoil equivalent energy: [2, 6] keVee and [4, 6] 
keVee- To convert the energy ranges defined in the ex- 
perimentally measured units into recoil energy for the 
appropriate target nucleus, we use q - 0.3 keVge/keV 
for Na recoils and q - 0.09 keVge/keV for I recoils, 
where Et-^vee - QErecoH and q is the quenching factor 
taken from [39|. The energy threshold of 2keVee is 
chosen based on the analysis thresholds and trigger ef- 
ficiencies achieved by previous experiments such as 
DAMA/LIBRA ITO! and NAIAD IHl. The higher 
4keVee threshold provides a comparison to a slightly 
more pessimistic scenario. The upper bound is 6 keVge- 
Above this, the rate for recoils, in a standard halo model, 
drops well below the considered background rates of 
0.5 - 5.0cpd/kg/keVee. In order to compare to DAMA 
data, we define the constant rate of events, A^o, which is 
the sum of the time-independent dark matter component 
(Ro) and the time-independent background. 

The sensitivity is calculated by fitting equation|3]to a 
simulated dataset containing no modulation signal in it 
(null hypothesis). The simulated data are created by ran- 
domly drawing events from a Poisson distribution with 
mean given by Nq. We report on several trials with A^o 



ranging from 0.5 to 5 cpd/kg/keVee- The value A^o = 0.5 
corresponds to an optimistic goal of achieving a better 
background rate compared to DAMA/LIBRA (No ~ 1). 
The highest rate (A^o = 5) is characteristic of the back- 
ground levels achieved by NAIAD ifTTl . To obtain a 
representative statistical uncertainty, the total number of 
thrown events corresponds to the average number ex- 
pected for a 250 kg exposure over two years, integrated 
over the experimental energy range. For a given value of 
A^o, the simulated experiment is performed for the two 
energy intervals listed above. 

A simple chi-square fit allows So and S „, to float 
while keeping the frequency, cj, and the phase, t^, fixed. 
Under the assumption of a null result, the 90% C.L. up- 
per limit to Sm is calculated from the error returned on 
the fit. The resulting sensitivities, in the form of 90% 
C.L. upper limits to the spin independent cross-section, 
are shown in Figure [3] For comparison, the calculated 
DAMA/LIBRA 90% and 99.7% allowed regions, with- 
out channeling Il40ll . are also displayed on the figures. 

The results demonstrate that low energy thresholds of 
~2keVee, or lower, are critical for obtaining adequate 
sensitivity to constrain the DAMA/LIBRA allowed re- 
gions. Figure [5] shows that with a 2keVee threshold, 
an experiment that has 5 times the DAMA event rate 
(and hence ~5 times greater background) provides sig- 
nificant constraints with two years of exposure. Con- 
sidering the sensitivity curves for the higher threshold 
scenario, it is clear that control of background becomes 
more important. With a 4keVee energy threshold, one 
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must improve upon the level of background achieved by 
DAMA/LIBRA in order to completely exclude the 3cr 
allowed regions at 90% C.L.. This result is understood 
because the expected rate of recoils fall off roughly ex- 
ponentially with increasing recoil energy. The feasibil- 
ity of lowering the energy threshold of the experiment 
will depend on achieving a balance between increased 
signal and background rates. 

Examining one set of sensitivity curves with vary- 
ing background and fixed energy threshold, it is ob- 
served that the sensitivity for a given WIMP mass scales 
roughly inversely with the square root of A^o- This is 
easily understood because the sensitivity of the annual 
modulation search scales with the uncertainty in the 
number of unmodulated events, which is mostly back- 
ground in this case. This demonstrates that the modula- 
tion analysis implicitly performs a background subtrac- 
tion. For similar reasons, the sensitivity of a Nal exper- 
iment will increase with the square root of the exposure 
time or mass. 

In this study, we have only considered sensitivities to 
spin-independent elastic scattering with a standard halo 
model. The ability to confirm or exclude a DAMA-like 
dark matter signal for alternative scenarios (e.g. spin- 
dependent scattering, inelastic models, etc.) is depen- 
dent on the kinematic constraints implied by the mod- 
els. Understanding the sensitivity reach for Nal under a 
variety of energy thresholds requires individual inves- 
tigation. The conclusions on relative sensitivity with 
scaling of background and exposure will remain largely 
the same. As a model-independent exercise, we eval- 
uate the statistical potential of a 2-year, 250 kg expo- 
sure to confirm a DAMA-like signal without consider- 
ing the nature of the possible underlying dark matter 
signal. Events are drawn from a Poisson distribution 
assuming the same rate, phase and period of modulat- 
ing events as observed in the DAMA/LIBRA 2-6 keVee 
energy region f9l. We find that the statistical power of 
the 250kg, two year, simulated dataset has a >5cr sig- 
nificance and hence could confirm a DAMA-like signal 
to that level of significance. 

4. Summary & Outlook 

In this paper, we have described a proposal for an an- 
nual modulation dark matter search in Antarctic ice us- 
ing Nal scintillation crystals. The ice provides a clean, 
quiet and stable environment, which is ideal for an an- 
nual modulation search. Furthermore, seasonal effects 
on background will be out of phase with experiments in 
the Northern Hemisphere, thus providing an indepen- 
dent check of the DAMA results. The sensitivity calcu- 



lation presented in this work demonstrates that verifica- 
tion or exclusion of the DAMA signal can be achieved 
with a two year deployment of a 250 kg experiment, as- 
suming a comparable energy threshold and similar con- 
trol of background as the DAMA experiments. 

Preliminary simulation studies suggest that back- 
ground from the Nal crystals and the materials used 
to instrument and encapsulate the crystals can be com- 
parable or lower than background rates achieved by 
DAMA. Radioactive background from crystals will re- 
main a dominant component. Future R&D will focus 
on growing low-background crystals, designing low- 
background encapsulation for pressure shielding in the 
Antarctic ice, and developing a deployment plan. 

A feasibility study on the remote operation of Nal 
crystals in South Pole ice is now underway. Two crys- 
tals, totaling 17 kg, were tested at the Boulby Under- 
ground Laboratory and later moved to the University of 
Wisconsin at Madison where they were encased in stain- 
less steel vessels. During the polar season 2010/2011, 
the units were attached to IceCube strings and success- 
fully deployed to a depth of 2450 m in the ice. Data are 
being gathered from the prototypes. Results from this 
demonstration will be available in an ensuing publica- 
tion. 
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